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1 Introduction 



The purpose of this paper is to construct a Rabi oscillation in four and five level systems 
which is a natural extension of two and three level systems, jl], [2] or [Hj. 

That is, we consider a model of an atom with n energy levels interacting with n — 1 
external (laser) fields and solve the Schrodinger equation to obtain the Rabi oscillation. 
How to solve it can be reduced to the calculation of an exponential matrix, which is very 
clear from the mathematical view point. However, to calculate it in the general case 
seems almost impossible. For n = 4 and 5 we exactly calculate it by force to get the 
Rabi oscillation. The unitary matrix corresponding to the Rabi oscillation may play an 
important role in constructing the generalized Walsh-Hadamard matrix (transformation) 
or other useful matrices in four or five level systems. 

This is a succession of the preceeding paper ^_8jJ in which we treated the special case with 
coupling constants being all equal. Our motivation of the work is to construct Quantum 
Computation based on many energy levels (of laser-cooled atoms in a cavity). See also 
|nj-[ni for related work. 

2 General Theory 

We consider an atom with n energy levels {(|A;), E^.) |0<A;<n — 1} which interacts with 
n — 1 external fields. As for the external fields we use laser fields with frequencies equal 
to energy differences of the atom and set = E/^ — Eq for 1 < A; < n — 1 . Moreover we 
assume 

El ~ Eq > E2 — El > ■ ■ ■ > En-l — En^2- 

We subject the atom to n — 1 laser fields having the frequencies Uk equal to the energy 
differences — Afc_i = E^ — E^-i- See the following picture : 
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The Hamiltonian that we use is written as 
H ^ 
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g^e'^^it+'t'l) 

g^Q~ii'^3t+4>3) 

An-2 (7„_ie^('^"-i*+'^' 



where gi, g2, - ■ ■ , Qn-i are coupling constants and 0i, 02, • " " > 4'n-i some phases. Here we 
have assumed the so-called rotating wave approximation (RWA) from the beginning. We 
would like to decompose the Hamiltonian. If we set 



gi(a)it+0i 



gj(t<Jlt+W2t+</>l+</>2) 



then it is easy to see 
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91 Ai 92 

92 As 93 



(2) 



9n-2 A„_2 9n-l 
9n-l A„_i 
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To solve the Schrodinger equation (we set h — 1 for simplicity) 

dt 

for U{t) being unitary we set U = VU (<S^ U = V^U), then it is not difficult to see 

i—U — 
dt 



{ Eoln + 



gi 

gi Ai - LJi g2 

g2 A2 - a;i - UJ2 ga 



gn-2 A„_2 — J2k=l 



M-1 



gn-1 A„_i — J2k=i 



Now if we set the resonance conditions 



Afe — Afc_i — cok for 1 < k < n — 1, 



then the equation above reduces to a simple one 



where 



C = C{gi,g2,---,gn-i] 



gi 
gi g2 
52 gs 



gn-2 gn-1 

gn-i 



The solution is easily obtained to be U{t) = e **^oe so the solution we are looking 
for is 

U{t) = V^U{t) = e-^*^oyte-itc_ 
Therefore we have only to calculate e"'**", which is however very hard task. 

In the following sections let us list our calculations for n = 3, 4, 5. 



3 Three Level Systems 



We only list the calculation in the case of three level systems. 



exp(— ztC) = exp < —it 



^ g, ^ ^ 

V 92 y J 



/ 3fCOS(v/gf+g|t)+g^ 



—i- 



•91 



. giSin(^g2+g2t) 



Sin(^g2+g|f) gigaCOSCy/gf+ggt)- 



■9152 



9^+9^ 



gig2COS(yg^+git)-gig2 . gaSinC^g^+g^f) g|COS(v/g;+g|t)+gf 

»i+f2 y^+gf 9^+92 y 



(10) 



see Appendix in ^ or p. As a result the solution we are looking for in the three level 
systems is 



U{t) 



-itEo 



I 



\ 



1 



-i{y)\t+4>\) 



-i(a;it+cJ2t+<^l+<^2) 



X 



/ g^COS(^/g^+gp)+g; 



9^+92 



— I 



giSin(^g^+g2t) 



9i92COS(^gf+g|t)-gig2 
9?+9l 



— Z 



giSinC^gj'+gif) gig2COS(^gf+g^t)-gig2 \ 



A/9f+9: 



9^+9^ 



^. g2Sin(^g2+g2t) 
V9?+9i 

.g2Sin(^g2+g2f) g2coS(^gf+git)+g^ 



C0S( JC/I + 



V9f+9: 



9!+9r 



111 
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4 Four Level Systems 

We continue our calculation in the case of four level systems. We want to calculate 



exp(— ztC) = exp 



-it 



^ ^ 

^1 5^2 

V yj 



Oil a\2 OlS 0.14 



«21 ^22 «23 024 
flSl O32 ^33 O34 
y a41 042 0^43 CI44 j 



(12) 



exactly. To make C a diagonal form let us calculate the characteristic equation 
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find solutions 


with 


good form 


(see j 



\'-{9i^9i^9'^y^'^9i9i. (13) 



formula). For that we set 

A = 9l^ (91 + 9z)\ B = gl + {g^ - g^f , 
then all solutions are given by 



Ai- - , A2- - , A3- - 
Let us rewrite these. It is easy to see 



A4 



• (14) 



_^fA + ^/B _ gigs 
M — = ^1 A2 — — ^ — , A3 



9i93 



A4 



-A. 



(15) 



2 ' ' A ' A 

Next we must look for orthonormal eigenvectors |Aj) (j = 1, 2, 3, 4) corresponding 



to eigenvalues above which are rather complicated. 



|Ai) 



9192X 
X92X 
(A^ - 9!)X 



g3(A^-g?) 



IA2) 



9293y 



A 

,2 \2^ 



V i9i - y)Y I 



lAs) 



X 
Y 



' -Xg2Y ^ 




^ 9192X 


A ^ 


, |A4)= 


-Xq2X 
(A^ - 9!)X 






\ X ^ 


by 


1 





(16) 



2 m-9l + 9i + 9l) + 9K9I + 9l - 9l)} 
1 

/2 {9li-9l + + 9i) + >?{9l + i7l - 91)} 
Therefore we obtain the orthogonal matrix 

iy=(|Ai),|A2),|A3),|A4)) G 0(4), 

which makes exp(— ztC) an easy form to calculate 



(17) 



exp(-2tC) = W 



-it\2 



-itXs 



V 



-it\4 



(1^ 



Let us calculate components of the matrix ()12p . From (fTTjl and (fTHjl a long calculation 
leads to 



an = 2gl [glX^cos{\t) + A^F^cos (^t)} , 



«i2 = «2i = —'^iXgl < (7iX^sin(At) + (^sy^sin 



ai3 = aai = 2g,g2 {(A^ - (7^)X2cos(At) + {gl - X^)Y^cos (^^t) 
aM = a4i = -2^^72 | ^i^3(^^-^i) x^sin(At) + Xigl - X')Yhm (^^t 

= 2gl |A2x2cos(At) + glY\os (^t) } , 
a,3 = = -2.,, I A(A^ - ^?)X^sin(At) + '-^^^^Y^^r. (^t) ' 
a24 = a42 = 2^2^3 {(A' - ^7D^'cos(At) + {gl - X')Y^— ^ ^^^^ - 
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cos 



V A 



ass = 2 |(A^ - glfX'cos{\t) + ^M^y^cos 



as4 = a4s = -2i\ [g,{\' - gffXhiniXt) + g,{gl - X'fYhin 



044 — 2 



A2 



X2cos(At) + {gl - X'fY^cos (^t) | . 



(19) 



As a result the solution we are looking for in the four level systems is 



/ 



U{t) = e 



-UEo 



-i(a;if+0i) 



X 



(20) 



/ \ 

Oil ai2 flis a'14 

^21 0^22 a2S a24 
OSl flS2 flSS 0,34: 

y 041 a42 a4s 044 J 

We obtained the explicit solution of Q under the resonance condition ©• This is the 
Rabi oscillation that we want, which is rather complicated. 



5 Five Level Systems 



We want to calculate 



exp(— ztC) = exp 



-it 



' ^1 ^ 

^71 (72 

(72 (7s 

(73 ^4 

V (74 



/ \ 

ail fll2 OlS fll4 fll5 

a2i a22 a2s ^24 ^^25 

^31 O32 ass ^34 ^ss 

O4I O42 0-43 0,44 '^45 

y 051 a52 ass ^54 c^ss y 



(21) 
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exactly. The characteristic equation is 



0= IAI5-CI 
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-93 


A 


-94 











-94 


A 



= A {A^ - {gj + gl + gl + gl)X' + {gfgl + 9^1 + 9hl)} ■ 
All solutions are given by 

Ai = 



(22) 



+ v/-4- 2v/S 



A2 



A + 2VB- yA-2VB 



, A3 = 0, 



_ ^Ja + 2VB-^A-2Vb _ ^Ja + 2^/B + ^JA-2^/B 

A4 — 5 A5 — , 



(23) 



where 



For 



A = gi + gl + gl + gl B = glgl + glgl + glgl 



orthonormal eigenvectors |Aj) (j = 1, 2, 3, 4, 5) corresponding to the eigenvalues 



|Ai) 
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^ 919293X ^ ^ 
M293X 
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(A^ - - gDg^X J y 

->^9i9293Y ^ 
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^{§-9l-9l)Y 
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9193 , 



(25) 
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where X and Y are given by 
X = 



2 M + 9l)A - 2B} A2 + 2{g! + gl - gl - gl)B 

1 



^ pB{{gl + gl)A- 2B + {gf + gl - gl - gl)\^} 
Therefore by the orthogonal matrix 



iy = (|Ai),|A2),|A3),|A4),|A5)) G 0(5), 



(26) 



we have a diagonal form 



exp(-2tC) = W 



-itX2 



-itXs 



-itX4 



V 



J 



Let us calculate components of the matrix (PT|) . A long calculation leads to 



(27) 



an = ^ + 2glglgl i^X\os{Xt) + X'Y\os (^^t) | , 
ai2 = = -2tXgig^gl |x2sin(At) + ^/BVhm (^^tj | ' 

ai3 = asi = + 2g,g,gl (^{X' - gl)X\os{Xt) + {B - X'gl)Y'cos (^^t) | , 

ai4 = 041 = -2iXgig2g3 |(A^ - gl - g^X'^siniXt) + \fB (y^- gl- gtj ^^sin [j^t 



ai5 = asi = c/1^25354 j + 2(A' - - (72')X'cos(At) + 2{B - X'gf - X'gi)Y'cos 



B 



022 = 2g^gl { X'X''cos{Xt) + BY'^cos { 



a23 = a32 = -22Af72^3' I (A' - ^D^'sin(At) + - gl^ Y\m | , 



a24 = a42 = 2g2g3 { X\X^ - gf - gl)X^cos{Xt) + B{^-gl-gi^ Y^cos l^t\ \ , 



025 = ^52 = -2iXg2g3g4 I (A^ ~ gl- gi)X'^sm{Xt) + Vb (^-^ - gl - gtj ^^sin 
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034 = a43 = -2tXgs I (A^ - (7?) (A' - gf - gl)X\m{Xt) + 



9i 



B 
A2 



^7? - 9l ) ^^sin ( 



a35 = = + 2g,g, - gl){X' - g\ - gl)X'cos{Xt) + - 9^ ) x 



a44 = 2A^(A^ -gf- g',yX'cos{Xt) + 2Bi--gi- gi Y'cos ' 



A2 



A 



i2 „2 „2\2 ta2„- / \ .\ , ,r5l ^ „2 „2A -i^2^- "^^^ 



a45 = a54 = -2zA^4 |(A' - gi - giyX'siYi{Xt) + V5 - - g^^ F^sin (^^t 
a55 = ^ + 2^4^ I (A^ - - glfX^cos{Xt) + A^ (| - ^7? - ^7^^) V^cos [^t) | 



(28) 



As a result the solution we are looking for in the five level systems is 

\ 

I \ 

ttii flia ai4 ais 



[/(t) = e 



fl21 ^22 ft23 024 ^25 

^31 C^32 ^33 '^34 ^35 

Q.4I (342 0^43 0-44 0.45 

y flSl O52 053 054 ass j 

We obtained the explicit solution of under the resonance condition This is the 
Rabi oscillation that we want, which is quite complicated. 



(29) 
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6 Discussion 



We have calculated e~***" with (jH} for n = 4 and 5 which is a generalization of the case 
n = 3 in ^ or inj. However, to calculate the general case is a very hard problem. One 
reason is as follows. In the general case the characteristic equation |A1 — C| becomes a 
general algebraic equation with degee n. By the Galois theory it is impossible to solve it 
in an algebraic manner. That is, we have no algebraic method to write down the solutions 

For the special case gi = g2 = ■ ■ ■ = Qn-i = (? we can calculate it easily, see [Hj. 

We also note that problems of this type occur elsewhere in Mathematics and Physics, 
see for example [10 . Thus our results may have broader applicability. 
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